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1 Introduction 

The fact that topological string amplitudes are closely related to certain 
holomorphic quantities in the physical superstring theories was known for 
some time. A practical incarnation of this relation is the recently discovered 
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[3J lU IS] gauge theory - matrix model connection. By now, this idea has been 
investigated using three rather independent approaches. 

In the original approach of [HI EJ E|, one starts from the open/closed 
string duality implied by a geometric transition, and computes gauge theory 
superpotentials using the fluxes of the dual geometry [H|-[IH|. The result is 
expressed in terms of the partition function of a certain closed topological 
string theory. On the open string side of the duality one relates the terms of 
the effective superpotential with the partition function of the holomorphic 
Chern-Simons theory (which describes the fields on the wrapped D-branes 
sourcing the geometry). Dijkgraaf and Vafa conjectured that the open and 
closed partition functions are identical. Since the computation of the open 
string partition function reduces to the computation of the partition function 
of a large N matrix model, this conjecture implies a certain relation between 
gauge theory effective superpotentials and matrix models This relation 
was further strengthened by the study of the underlying geometry of the 
matrix model and of the gauge theory [U 

In the second approach [20], the effective glueball superpotential of an 
Af = 1 theory with adjoint matter was evaluated using superspace techniques. 
It was found that only zero momentum planar diagrams contribute to this 
superpotential, thus validating the original conjecture. 

In the third approach [HJ the generalized Konishi anomalies of the 
field theory were used to obtain relations between the generators of the chiral 
ring of the theory. These relations, which under certain identifications can be 
reproduced from a matrix model, can then be used to construct the effective 
superpotential. 

Perhaps the most immediate extension of the matrix model-gauge theory 
relation is to theories with fields transforming in matter representations of 
the gauge group j2S]-jlZl- For theories with fields transforming in the funda- 
mental representation of the gauge group it was suggested that the addition 
to the original DV proposal involves matrix model diagrams with a single 
boundary. For arbitrary generalized Yukawa couplings and a simple super- 
potential for the adjoint field this proposal was proved in Furthermore, 
it was shown in |7j that the matrix model fully captures the holomorphic 
physics of theories with Nf < N c , regardless of the complexity of the tree 
level superpotential for the adjoint field. 

Extending the correspondence to gauge theories with baryons turned out 
to be somewhat more challenging. In particular, baryons only exist for the- 
ories with certain relations between Nf and N c ; therefore, taking the large 
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iV c limit to restrict to planar diagrams is not possible. Moreover, the no- 
tion of boundary in diagrams with baryons is not well-defined unless certain 
manipulations are performed. 

In [32] an extension of the DV proposal to theories with baryons was 
formulated and it was shown that for an SU(N C ) theory with Nf = N c 
quarks this proposal reproduces exactly the known gauge theory physics. 1 



The goal of this paper is to extend the correspondence to supersymmetric 
QCD with Nf = N c + 1 flavors. As it is well known [2], this theory has Nf 
baryons, Nj mesons and a dynamically generated superpotential 

W = J ± z - 1 (B l Mi&-detM) . (1) 

In order to relate this theory with the corresponding matrix model, we deform 
it with the appropriate sources, and integrate out the mesons and baryons. 
We then compare the result with the one given by matrix model planar one- 
generalized) boundary and find perfect agreement. 

Our interest in theories with baryons has several reasons. The first is that 
neither of the three routes by which the original matrix model-gauge theory 
relation was reached seems easily extendable to these theories. 2 Second, 
in theories with only chiral flavors, baryons are the only objects one can 
construct. If one is to extend the matrix model-gauge theory relation to 
such theories, understanding the role and the correct treatment of baryons 
is crucial. 

Last but not least, the superpotentials of theories similar to SQCD with 
Nf = N c + 1 flavors 3 are used as starting points for the construction of low 
energy effective superpotentials in many theories where symmetries do not 
determine these superpotentials directly. It is therefore important to have a 
direct method of computing them. 

1 Related work has appeared in |3Uj . where a perturbative field theory computation in 
the spirit of |20j was used to recover the terms linear in baryon sources in the gauge theory 
effective superpotential. 

2 In particular, as we will see in the last section of this paper, understanding the baryons 
at Nf = N c + 1 in geometric transitions is quite difficult. 

3 These are the so-called s-confining theories. They have a description in terms of gauge- 
invariant composites everywhere on the moduli space, and the effective superpotential for 
the confined degrees of freedom is not singular at the origin of the moduli space. (For 
discussions on s-confining theories see 0^1 EDI an d references therein). 
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One of the ways in which the validity of the Nf = N c + 1 superpotential 
is usually tested is by obtaining the correct Nf = N c superpotential in the 
absence of baryonic sources. However, this superpotential contains much 
more information, which can only be captured by turning on all the baryon 
sources. The fact that the matrix model reproduces the rather involved 
effective superpotential obtained with all the baryon sources turned on is 
a very powerful confirmation of the validity of the extension of the matrix 
model-gauge theory relation to baryons. 

Integrating out all the fields in both the gauge theory and the matrix 
model is rather complicated, and often results in rather unedifying expres- 
sions involving roots of large degree polynomials. Fortunately, there exists a 
procedure j2Zj which allows us to compare gauge theory and matrix model 
results in theories containing only mesons. Thus, to compare the gauge the- 
ory and the matrix model results it is enough to integrate out only two flavor 
fields on the gauge theory side (sections El and EJ), and to relate the result to 
the matrix model free energy obtained by treating the N c — 1 massless flavors 
as background fields (section 0). 

Except for SU(2), where computing the values of the effective superpo- 
tential at its critical points is not too difficult (sections 2 and 3), we will only 
be comparing matrix and gauge theory results at Nf = N c — 1. In section 
18.21 we will derive these values, using gauge theory techniques, and discuss 
a method for computing the full matrix model free energy. As an example 
we apply this method to the case of an SU(2) gauge theory and recover all 
expected results, already discussed in sections l3~2l and 141 

Before proceeding, let us remark that in our case there is no distinction 
between the unitary matrix model and the hermitian one. This is due to the 
fact that we will be interested in theories containing only fields transforming 
in the fundamental representation. Since these fields are not constructed out 
of generators of the gauge group, they are the same both in the SU (N c ) and 
in the U(N C ) theories. Thus, the matrix integral is the same for both gauge 
groups. 

To fix the notation, Latin indices from the beginning of the alphabet 
(a, b, c) are SU(N C ) indices; Latin indices from the middle of the alphabet 
k) are SU(Nf = N c + 1) indices; Greek indices from the beginning of the 
alphabet (a, j3, 7) are SU{2) indices corresponding to two flavor fields which 
are singled out. They take the values N c and N c + 1. Hatted Latin indices 
from the middle of the alphabet are SU(Nf = N c — 1) indices, correspond- 
ing to the flavor symmetry unbroken by the quark masses (but nevertheless 
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broken by the presence of baryon operators). 

2 Review of the Dijkgraaf-Vafa proposal for flavors 

In a series of papers [HI IH E|, Dijkgraaf and Vafa proposed a perturbative 
method for computing the effective glueball superpotential of certain Af = 1 
theories with fields transforming in the adjoint and bifundamental represen- 
tations of the gauge group. According to this proposal, the planar diagram 
contribution to the free energy of a certain matrix model yields the effective 
superpotential of the corresponding M = 1 gauge theory. In 't Hooft's double 
line notations these diagrams have the topology of a sphere. 

When fields transforming in the fundamental representation of the gauge 
group (quarks) are present one must also include the free energy arising from 
planar diagrams with one boundary (diagrams with the topology of a disk) 
[231I2S1- More explicitly, the gauge theory effective superpotential is given 
by 

W cS (S, A) = N C S(1 - In — ) + N c -£± + N f T x=1 , (2) 

where the first two terms are also present in a theory with only adjoints, and 
the third term is the contribution of the flavors. 

If baryonic sources are also added, the diagrams that can be constructed 
become more complicated. However, only the planar diagrams with as many 
index loops as a diagram with one boundary contribute to the effective su- 
perpotential Since the number of colors and the number of flavors are 
related, it is not possible to select the relevant diagrams by taking the limit 
in which the number of colors is large. Thus, the planar diagrams with 
baryon sources have to be selected by hand. In |32] ^ was shown that planar 
baryonic diagrams for SU(N C ) theories with N c flavors reproduce the known 
physics. 

The nonlinearities introduced by the baryonic operators make the compu- 
tation of the matrix model partition function challenging. Fortunately, when 
the tree level superpotential can be expressed in terms of mesons, the matrix 
model and gauge theory can be related more directly [27|. Thus, adding 
in the matrix path integral a constraint which identifies the matrix model 
quark bilinears with the gauge theory mesons allows one to compute directly 
the gauge theory superpotential with the corresponding mesons integrated 
in. This proposal was proved using the geometric construction of the matrix 
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model [TBI, using the symmetries of the gauge theory [22], and by explicitly 
integrating in quarks [¥T] . 

Thus, the free energy which gives the superpotential of a theory with 
both massive and massless flavors is given by 

e T = I DQ DQ5(Q~Ql - Mhe~ w ^ QM (3) 

J v ' " 11 planar+l boundary v ' 

where Q a are the massive quarks and Q% are the massless ones. 

This expression also allows us to compute the free energy in the presence 
of baryons, as long as Nf — N c + 1. Indeed, choosing two massive quarks, 
it is possible to sum up all Feynman diagrams involving them and obtain a 
result which only depends on quark bilinears. 4 At this stage the 5-function 
constraint can be easily enforced and we are left with computing the integral 
of the constraint. 

If the dimension of the matrices M c is larger than Nf (which is always 
the case in large M c limit) this integral is 



DQ i DQ j 5(Q CL Q^ - M?) = e M cM&ct m / k 2N f )-n f \n(& c t m / k 2N f )-n s M c \nM c ^ 

(4) 

where A is a cutoff. This result however contains both leading and subleading 
terms in Nf/M c . In particular the logarithm in the second term in the 
exponent is proportional to the number of flavors, Nf (the determinant is 
taken over flavor indices, and therefore is of order M N f), and therefore this 
term is of order Nf. Hence, it is generated by a multi-boundary diagram 
with insufficient gauge index loops, and should not be included in the free 
energy. 

Identifying the matrix model 't Hooft coupling with the gauge theory 
glueball superfield and taking into account the clarifications above, the con- 
tribution of the planar and 1-boundary diagrams to the above integral be- 
comes: 

DQ t D&5(Q?Ql ~ = e SH*etM/A™f)-N f SHS/AS) (5) 

This equation will be one of the important ingredients in our comparison of 
matrix model and gauge theory results. 



4 This little miracle happens only for Nf — N c or Nf — N c +1. If Nf > N c + 2 one needs 
to choose more than two massive flavors, which makes the matrix integral non-Gaussian 
and rather hard to compute. 
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3 Integrating- Out All Flavors 



In this section we compute the gauge theory effective superpotential at its 
critical points. The unbroken symmetries determine its value up to an un- 
known function of one variable. By requiring consistency with the high en- 
ergy theory, we construct a differential equation for this function. We first 
solve it for the special case of an U(2) theory with three flavors and then 
turn to the general case. 

3.1 Symmetries and Consistency Constraints 

We start with a tree level superpotential with mass terms and baryon source 
terms for all flavors: 

W tree = m)Q^ + b % B* + VBj . (6) 

The total superpotential is the sum of this tree level superpotential and 
of the dynamically generated superpotential 

W *»» = "75^=1 ( det M - BM ®) ■ ( ? ) 



The quantum numbers of the sources are 
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Using these quantum numbers, we can determine the form of the allowed 
superpotential terms after integrating out all flavors. 

All superpotential terms are functions of A 27Vf ~ 3 , b, b, and m. To con- 
struct invariants under the non-abelian flavor symmetries, the only allowed 



7 



building blocks are bim^W and detm. The U(1)r invariant combination 

of these is (6 i m}^') JV /- 1 /(det m) 2 . Its U{1) A charge is (-2N f )(N f - 1) - 
2(-2)N f = -2Nj + 67V/. Therefore the combination 



(&im#) 



JVf-l 



(detm) 2 1 1 K ' 

is invariant under all symmetries, and is dimensionless as well. 

The existence of the gluino condensate implies that in the absence of bary- 
omc sources, o = 6 = 0, the superpotential is (N f - l)[(detm)A 2N f- 3 ] 1 ^ N f- 1 \ 
Therefore the possible form of the superpotential in the presence of baryon 
source terms is 

W eff = (N f - l)[(detm)A^-Y^/ ( (A^Y^ , (9) 

where f(x) is a function we want to determine. 

In the limit of infinite mass parameter m this theory reduces to a pure 
N — 1 gauge theory. In this case we know that there are N c vacua and the 
values of the superpotential at the critical points differ by roots of unity of 
order A^ c . In this limit, the argument of the function / in the equation above 
vanishes, while its coefficient can be identified with the dynamical scale of 
the resulting theory. Thus, in order to recover the expected gauge theory 
results, we must impose the boundary condition 

/(0) = u% c , k — 0,...,N c — l with u Nc = l . (10) 

The expectation values of the moduli are obtained by differentiating this 
effective superpotential with respect to the sources. 5 

Mi = ^: = [(detm)A 2 ^- 3 ]W 1) x 

(N f - l)xf(x) { ^ Nf ^ lV ~ 2(-- 1 r 
= ^^ = (AT / -l) 2 [(detm)A 2JV /- 3 ] 1 /( Ar /- 1 )x/ / (x) m ^ 




(bmb) 



B 3 = ^Z^ = (N f -l) 2 [(detm)A 2N f~ 3 ] 1 ^ N f-Vxf\x)-^ 



(bmb) 



3 We use the simplified notation (biinW) = (bmb) 
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Therefore, we find 



B*Ml = (N f - l) 2 [(detm)A 2N f- 3 f (N f-^xf'{x)x 

( (m-')i/(x) + (N, - !)*/<(*) ( {N <- ^ - 2( m -^ 



(bmb) J J (bmb) 

= (N f - l) 2 [(detm)A 2N ^ 3 ] 2 ^ N f^xf'{x)-^- x 

(bmb) 

(f(x) + (N f -l)(N f -3)xf'(x)). (14) 

One of the equations of motion derived from W tre e + Wd yn (0 0) imposes 
the following relation: 

ir\i; .\ 2V - . (15) 

Therefore, 

(N f - l) 2 x- 1/{N i- l) xf(x) \f{x) + (N f - l)(N f - 3)xf'(x)\ = -1. (16) 

This equation is special for Nf = 3, as the secont term vanishes. We 
begin in the following subsection by analyzing this case, and defer the general 
discussion to section 8. We then give the matrix model description of the 
SU (2) theory and compare the results. 

3.2 SU(2) with three flavors: Field Theory 

For Nf = 3 the superpotential © is 

M^2 l(d e t „,AT'v(gg). (IT) 

Let us notice that the argument of f(x) does not depend on A. This is easy 
to understand. Because of the Lie algebra identification SU(2) ~ Sp(l), the 
baryons can be interpreted as mesons in the Sp(l) theory. The mass matrix 
is 

/ e^ k bk m) \ 



-m) e ijk b k 



(18) 



-ijh 

and its Pfaffian can be perturbatively expanded around detm in inverse 
powers of m. The function f(x) must be precisely this expansion. Therefore 
it is a polynomial of y/x. 
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Indeed, equation (|TB|) reduces to 

Ax 1 l 2 f\x)f{x) = -1. 
which implies that f(x) is given by 

f{x) = ±{c-x^y\ 



(19) 



(20) 



The integration constant C is fixed to unity by the boundary conditions ()10|). 
Then, the effective superpotential becomes 



W = ±2[(detm)A 3 ] 1/2 ( 1 - 

1 1 detm 



1/2 



±2 



(detm) — (bmb) 



1/2 



A 3/2 



(21) 

The combination in the square bracket is precisely the Pfaffian of the mass 
matrix including the baryon source terms. In the next section we will recover 
this result from matrix model computations. 



4 SU{2) with three flavors; Matrix Model 

In this section we describe SU(2) supersymmetric QCD with 3 flavors using 
the matrix model. Since the baryon operators in this theory are bilinear in 
quarks, the matrix model free energy can be computed directly. We will find 
that, after integrating out the glueball superfield, the effective superpotential 
agrees with the field theory result given in Eq. ()21j) . 6 

As briefly stated in the previous section, for an SU (2) theory with three 
flavors the tree level superpotential is 

W tree = miQIQi + b^ k e ah Q a 3 Q\ + V e l]k e ab QiQ k b . (22) 

To compute the partition function it is useful to rewrite this expression as 

W tree = ^Q T K u(2) Q (23) 



where 



Kui ^{ m!e« tC«^J WJ 1 (24) 



6 For larger N c the matrix model is sufficiently complicated to render challenging the 
direct recovery of the field theory results. We will return to these questions in section [HI 
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The 1-boundary free energy is given by the logarithm of the determinant 
of if j/(2). This can be easily computed and it gives: 

detK u{2 ) = (detm - (bmb)f X2 (25) 

where in the exponent the first factor of 2 is due to the fact that we integrated 
over two types of fields, Q and Q, while the second factor of 2 represents the 
number of colors. 

In principle one should worry about isolating the planar diagram con- 
tribution to the free energy. Fortunately, for SU(2), all the diagrams are 
planar. 7 

Combining this with the Veneziano-Yankielowicz term yields the effective 
superpotential: 

W eff = N C S (l - In — \ - Sin (det m - (bmb)) (26) 

To compare with the field theory result we must integrate out S. This 
gives 

W eg = ±2 (det m - {bmb)) 1/2 A 3/2 (27) 

which precisely matches the field theory result. 

Perhaps this agreement should not appear surprising, since for a U(2) 
gauge group the mesons and baryons have similar structure. However, the 
computations which led to the two results are substantially different; this 
seems to imply that the agreement is somewhat nontrivial. Another point 
worth emphasizing is that all matrix model diagrams contributed to the 
effective superpotential. The origin of this fairly surprising fact is again the 
bilinearity of the baryons. This will not happen in the general case to which 
we return in section [7| 



5 Integrating-Out Two Flavors - The Elegant Way 

As discussed in section our goal is to match the gauge theory effective 
superpotential after integrating out two quarks with the matrix model pre- 
dictions. Let us therefore begin with the appropriate computation on the 

7 This is due to the fact that both 6^ as well as e a b are invariant tensors. The non- 
planarity can in principle arise due to insertions of a baryonic operator in the Feynman 
diagram, but the antisymmetry of e a {, can be used to transform it into a planar one. 
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gauge theory side. There are two ways to achieve our goal. In this section, 
using symmetry arguments, we constrain the form of the effective potential 
after integrating out two quarks and then derive certain constraints on the 
unknown functions. Solving these constraints leads to our result. In the next 
section we rederive the same result by directly integrating out the appropriate 
fields. 

Since we only give mass to two of the flavors, the tree level superpotential 

is 

W = m a p Q a QP + b t B l + VBj , (28) 

where a, f3 — N c , N c + 1. Hereafter we distinguish the indices of the massive 
and massless flavors: a, (3 — N c , N c + 1, and i, j — 1, • • • , N c — 1 respectively. 
The superpotential has a tree level part 

W = m a p M^ + b a B a + b i B i + b p Bp + b j B j , (29) 

and a non-perturbatively generated part 

(S'M/Sj + B a MlB s + B'M^Bp + BTM^Bp - det M) . (30) 
5.1 Preliminaries 

Our goal is to find the effective superpotential after integrating out the two 
massive flavors. This superpotential is a function of b a , b^, Jr, V, m and M~ . 
In the absence of baryonic source terms, it is easy to find the solution 

B a = B i = Bp = B s = , (31) 
Mp = = , (32) 

m_ 

detM 



(m' 1 )g(detm)A 2iVc - 1 
Mp = ;/3 \ ■ J- , (33) 



where, (m" 1 ) is the inverse of the two- by-two mass matrix, and M is the me- 
son matrix constructed out of the remaining flavors. The resulting effective 
superpotential is 

w= (detm)A^ 
det M 

which is the expected Affleck-Dine-Seiberg superpotential. 
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In the general case the quantum numbers under the SU(Nf — 2)q x 
SU(2)q x U(1)q C SU(Nf) flavor symmetry and its counterpart for Q, force 
the superpotential to take the form 

(detm)A 2A ^-i . 

W = K f v /(x,y), (35) 

det M 

where the invariants x and y are 

M)(detM) 2 

(detm) 2 A 2 ^" 1 ' 1 J 

(bM~ l b) detM 
* = (detm) • (37) 

with (bmb) = 6 a m^6 /3 and (bM~ l b) = ft^M -1 )^. We require that the 
superpotential be regular in the limit of no baryon sources and also at weak 
coupling A — > 0. Therefore, the function f(x,y) can be at most linear in x, 
and hence 

(detm)A 2JVc " 1 . . (6m6) det M , . . , . 

W = - g(y) + - - 1 h(y) . 38 

detM yyU ' detm yy> 1 ; 

In order to obtain the explicit forms oig(y) and h(y) it is useful to consider 
several limiting cases. 

5.2 bt = V = with i,j = 1, • --,iV c - 1 

In this case, the SC/ (Aj — 2)q x S77 (Afy — 2)q symmetry is unbroken. Hence, 

Mp = Mf = 0, (39) 
B l = B 3 = 0. (40) 

The equations of motion are 

B a B(s - (M _1 )£(det M) + A 2JVc_1 m^ = 0, (41) 
5" = (M" 1 )^, (42) 

= ^(M" 1 )^. (43) 

where {M~ l ) a ^ is defined only in the two-by-two block. Substituting the 
solutions from the last two equations into the first one, we find 

[{M-^fb^iM- 1 )}] - (M _1 )^(det M) + A 27Vc_1 m^ = 0. (44) 
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This is an equation for two-by-two matrices and hence there are four un- 
knowns. On symmetry grounds we take the following ansatz: 



(M )p = am* + (3(mb) a (bm) p (45) 

where a, f3 are function of the invariants. Apparently this system is over- 
constrained, as there are four equations for two unknowns. However, a solu- 
tion exists and is given by 

a = , -( de ^)( d ^) , (46 ) 

-(6m6)(detM) 2 + (detm^A 2 ^" 1 

P = . . (47) 

(det m)((bmb) (det M) 2 - (det m) 2 A 2 ^- 1 )A 3 ~ 27V / 

Using this solution, the superpotential is given by 

(detm)A 2jVe - 1 _ (det M)(bmb) 
eff ~ (det M) (detm) ' 

which is precisely what we expected from the symmetry considerations, ex- 
cept that we now determined the coefficient —1 for the second term. This 
determines the boundary condition h(0) — — 1. 

5.3 b a = P = with a, (3 = N c , N c + 1 

The next simple case is b a = \f = 0, when the only parameters in the effective 
superpotential are mf, bi, V , Mi. Hence there are no doublet breaking 
parameters of SU(2)q x SU(2)q. This immediately gives 

Mf = Mf = B a = Bp = 0. (49) 

The equations of motion can be easily solved, 

B s = -A 2Nc -\(M-% (50) 
B- t = -K 2N *-\M~%V, (51) 

M p -{m ), (det ^ } , (52) 
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where M is the meson matrix for the remaining Nf — 1 flavors. Substituting 
the solutions to the superpotential, we find 

= (detm)A^- 1 _ j^-i^-i^a (53) 
detM y ' 

and therefore g(y) = 1 — y. The only remaining function to be determined 
is h(y). 



5.4 General Case 

Putting together what we have learned so far, the superpotential is 

(detm)A 2JVc ~ 1 k2N Ul _- lZ . (bmb) det M 7 . . 

W eff = - '— A 2Nc 1 (bM 1 b) + - / h(y), (54) 

B detM 1 ^ detm Kyh 1 ; 

with h(0) = — 1 and 

(bM- l b)detM 

y = — u i \ — • 55 

(det m) 

From this superpotential we can obtain the vacuum expectation values of the 
M@ mesons and of the baryons: 

I, = dW eff _ (m-^(detm)A 2 ^ 1 (bjf) det M 

detM detm W> 

(bmb)detM (bmb) detM _ ± * 

%) n yh{y){m )£, (56) 



B* 



dW eff 








— (m~ 


J a 


dW eff 








dW eff 




db? 





detm detm 



, - - -A^-\bM-i) j+ ^ detM yh>(y)^i, (57) 

db 5 V ; ' detm y yy> (bM-tb) V ; 

<-'"<u (Ml detM 
Bp ~ del^ %) ' (58) 

The one piece of information we cannot obtain from this superpotential is 
the vacuum expectation value of the off-diagonal mesons. By symmetry 
considerations it must be of the form 

Ml = a(x,y)bJP, (59) 
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and similarly for M~ . To determine the unknown functions a(x, y) and h(y) 
one must use the equations of motion. We start with 



dW _ M?B P + M' a B. 



dB a 



J\2N c -l 
, 2N c -li 



'+b a 

(bmb) (det M) 2 



= j^-i{ h ^ N ^y) + b ^etmy J yI/i IJ),I(!J) 

(bmb) det M , 



+a 



-A 



2N C -1, 



(6M" 1 6)6 a + 



det m 



-yh(y)b a 



+ b 



This leads to the differential equation 



1 + h(y) + xyh'(y)h(y) + a—-^[-y + xyh'(y)\ = 0. 

det M 



(60) 



(61) 



Another useful equation is 

dW M?B B + MfBj 







dB i 



(bmb) det M 



\2N C 



a 



+ 



det m 

{bmb) det M 



h(y) - A 



27V C -1 



det m (bM 

This leads to another differential equation 

det m 



^tyh'iy^+k. (62) 



a- 



detM 



%) + h'(y) = 0. 



(63) 



Solving for a from the second equation and substituting it into the first 
one, we obtain 



1 + h(y) - xyh'(y)h(y) + (y - xyh'(y)) 



h'(y) 
%) 



0. 



(64) 



Because this equation has to hold for any x, it gives two equations for h(y), 

(65) 
(66) 



i + % )+y £M = o, 

h'(v) 

-yh>(y)h(y) + yh>( y )^L = 0. 



Ky) 
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It is non-trivial that two different non-linear differential equations have a 
consistent solution. The first equation gives 



dh 



h 2 + h 
1 + %) 



dy 

5 

y 



(67) 



and hence 

1 4- hU,\ 

log \y\ + const. (68) 

Together with the boundary condition h(0) = —1, this leads to the solution 

1 



%) 



On the other hand, the second equation gives 

dh 



and hence 



%) 



(69) 



(70) 



(71) 



Both equations give the same solution, which confirms our result. 8 There- 
fore, the effective superpotential after integrating out two quarks is 



W eff = A 



2JV C -1 



det m 
detM 



(bM- l ~b) 



(bmb) det M 



detm- {bM~ l b) detM 



• (74) 



6 Integrating-Out Flavors - The Laborious Way 



In this section we will recover the results of the previous section using a 
different method: instead of using symmetries to constrain the final form of 



3 We can also determine a(x,y): 



a(x,y) 



det to 



det M l-y 



and hence 



Ml 



b a V detM 



M9= ¥" detM 
3 1 — y det to 



1 — y det to 

These expressions are necessary for integrating the two flavors back in. 



(72) 



(73) 
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the effective superpotential, we will just directly solve the classical equations 
of motion and then evaluate the initial superpotential at these values of the 
fields. Starting from flU EOJ) 

W eff = M^rrip + bjB i + UB^ + b a B a + b a B a (75) 
+ -^-j [B^jBs + B l M«B a + B a M 3 a B 3 + B a M^B„ - det M . 

the equations of motion are: 

11, ll 1 Mf ■ ////.\ 2V ! = (76) 
BiM' + bjA 2 ^- 1 = (77) 

where i, j = l...N c + 1, M- is the cofactor, and only mP a ^ 0. We split the 
(N c + 1) x (JV C + 1) matrix M\ into a 2 x 2 block M% , and a (N c - 1) x (JV C - 1) 
block M~. The off diagonal blocks are M| and Mf respectively. 

Multiplying (J7|) by M{ and using the fact that M?M} = det M5\ we find 
after a few straightforward steps: 

Mfm^ = biB p (78) 

JW>£ = b% (79) 

M 7 Q m^ = 6 7 5^ + 8% det MA" 12 ^" 1 ' (80) 

M°m} = b a B p + <f£ det MA^ 2N ^ (81) 

Equations (|8U|) and (J81)) give Bpb a ml = B^bam^, which implies 

£7 = B ml ^A 2 ! 2 ^- 1 ) (82) 
5 a = B ml b^ 2N ^ , (83) 

where B is a parameter. 

We will first express all the expectation values in terms of B, and then 
use some of the remaining equations of motion to relate B and det M. The 
mesons are given by: 

Mp = b p b a BK 2{2N ^ l) + (m-% det MA- (2Arc - 1} (84) 
Ml = b~b l Bk 2{2Nc - 1] (85) 
M? = hb a BK 2{2Nc - l) (86) 
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Combining these equations with equation (J77|) one finds the baryons B l : 



B l 



— (M~ 1 )l6 J (l + X.B 2 )A 2JVc-1 
-b s (M- 1 ){(l + XB 2 )A 2N '>- 1 



where 



Substituting everything back into W e g we find 



(87) 
(88) 

(89) 



+ 



-K 2Nc - l (bM- l ~b){B 2 X + I) 



(B 2 X + 1) det M + BX(B 2 X + 3) 



y\2JVc-l 



(90) 



where as before (6M _1 6) is a shorthand for 6?(M _1 )j^. 

The next step in our evaluation is to find the relation between det M and 
B. Using the block decomposition of the meson matrix we outlined in the 
beginning, it is not hard to find that det M can be expressed as: 



det M = (det M"|) det(Mf - M^(M -1 )jMj 



(91) 



After expressing all its components in terms of B, one can easily compute 
the determinant of the 2x2 matrix to be: 



det (Mf-M^M" 1 )^ 



det M 2 + X \B - A 2 ^ 2N ^B 2 (bM^b) 



detM 



A 2(27v c -i) det 



ni 



(92) 

We should note that if k flavors were integrated out, the numerator on the 
right-hand-side of the equation above should be replaced with det M k + 
X(B - B 2 bM~ 1 b) detM**- 1 . Thus, the first equation relating B and detM 
is: 



(det Mf 



(detM) (det m)A 



2(2JV C -1) 



(det M) 



(det M)X{B - A 2 ^- 1) B 2 {bM~ 1 b)) 

(93) 

To find the other relations between B and det M we use the equation of 
motion: 



B a B p + m%A 



a a 2N C -1 



ii 1 



9 detM ,^ adet 
(detM/) 



Mf - M^M-^JMI 



(94) 
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Multiplying this equation by (Mf - M^M^lM?) and using the fact that 



Mf 



ddet Af 



det M <5f we obtain after a few steps: 



det M = -BX - - + A 2 ( 2iVc - 1 )(fi 2 X + lX&Ar^) . 



(95) 



Again this relation is independent of the number of flavors integrated out. 
One can also evaluate by hand the cofactors in (J94j) , sum them with mr a , and 
obtain 

j j. A2(2JV C -1) y 

det M = , (2 + B 2 X) - -(B - A 2 ^ 2N ^ B 2 {bM~ l b)) (96) 

2detM 1 J 2 K 1 

The equations ESI EEj) have a unique solution 

detm A 2 ^" 1 



detM : 

B -■ 
which gives 



-X 

detM 

A -2(2iV e -l) det M 

(bM^b) det M — det m 



B-A 2( - 2N ^B 2 {bM- i: b) 



det m 
detM 



(bM^b) 



(bmb) det M 



detm - {bM- l b) detM 



(97) 



(98) 



We have thus recovered the effective superpotential (|73jl constructed in sec- 
tion EJ We now turn to the matrix model analysis of the theory and recover 
the same results. 



7 SU(N C ) with N c + 1 flavors; Matrix Model 

The tree level superpotential of the theory under consideration was described 
in section 13.11 Since the goal is summing all diagrams containing two flavor 
fields, it is useful to rewrite it in the following form: 

W tree = m p a Q a p Ql + e aP b a Q a p V a + e a(3 b a Q^ a V a 

+ b^QlQ\Vl h + Ve af3 Q<* a Q$Vf b 

where a and (3 take the values N c and N c + 1, and 

\f — JV c ,JV c +l,si,...,!iv c -i f^a\ r^ a N c -i 

v a — t t aai ... ajVc _ 1 ^i 1 • • • Vj^-i 

ii,...,iiv c _i = l,...,JV c -l a 1 ,...,a Nc _ 1 = l,...N c (99) 
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and similarly for V a . Also, 

T/» _ <r N c ,JV c +l,i,ii,...,ij Vc _ 2 , /Q ai f) a N c -2 

v b 1 b 2 ~ c t Ol02ai---ajv c -2V« 1 ■ ■ • Vj JVc _ 2 

?i, . . • ,?at c _2 = 1, • • • ,N C - 1 , 61, 6 2 , ai, • • • ,ojv c -i = 1 • . .iV c (100) 

and similarly V i ab . 

For latter convenience let us point out that: 

Q\V a = Q\V a = (101) 

where % and j take values only from 1 to N c — 1. There is no constraint of 
this sort for Q^V^ b , etc. However, one can see that 

KW 6 = and V a vHl = (V) i = l,...,N c -l . (102) 

To systematically compute the integral it is useful to write the tree level 
superpotential as a quadratic form. This is easily done by introducing 

Q ~\QU v -\^A\e^Pv a ) (103) 

and 

(e aP ®b iV : ia2 mi® 5JI ) K T = Ke=(° M (104) 

Then, the tree level superpotential — the matrix model potential — can be 
written as: 

W tree = ^Q T KQ+Q T -V 

= ^(Q + K^vf K(Q + K- l V) -^K^V . (105) 

Therefore, the partition function is 

Z= f DQ i DQ j 5(QfQi~M^ vTK ' lv -^ ladetK . (106) 

J J K 1 01 % ' planar+l boundary v ' 

The exponent of the integrand can be easily analyzed; fairly standard 
matrix manipulations lead to: 



det if 



det c (5 b a det m + bftV^Vf) f , (107) 
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where det c denotes a determinant over the color indices a, b, while introduc- 
ing the notation V = hV^ and similarly for V, the inverse of K is given 
by: 

K - 1= (-e®V(t c detm + VV)-' m -i g, J!Z_)-i \ 
V m- 1 ®(l c + ^)- 1 -e^Vitcdetm + VV)- 1 ) 

Let us now analyze in some detail the combination (VV + l c detm) _1 . 
Equation (jlUlj) implies that we need to compute only the terms proportional 
to the identity matrix. The other terms will vanish upon contracting with 
V. It is not hard to see that 

(VV) b a = -bM-'b det M8 b a + Xi(M)Q«Qi, (109) 

where we have already used the 5-function constraint from the path integral 
to replace Q^Qi with M{. This in turn implies that 

(VV + det m)" 1 b a = ° + Y?Q?Q J 6 . (110) 

J det to- (6M- 1 6)detM J 

The precise value of Y is irrelevant, since the last term always cancels due 
to contractions with V a or V b 
Thus 

V*K-i V = ^»K= . (ill) 

det m - (bM~ l b) det M det m - (bM~ l b) det M 

Combining all pieces together we find that the gauge theory effective 
superpotential is given by: 

-Weff = S{l-In^) + Sln^^ (112) 
+ M)detif , - In det ( det m5 b + kPVLVf 



detm - (bM~ l b) detM 

The unit coefficient in front of the Veneziano-Yankielowicz term arises as the 
difference between the number of gauge theory colors iV c and the number 
of massless flavor fields Nf. Before we proceed, let us point out that the 
last term in the equation above has an implicit dependence on the glueball 
superfield. Indeed, as the determinant is taken over the matrix model color 
indices, the argument of the logarithm is of the order m . Exposing the 
part arising from the relevant planar diagrams is potentially complicated; we 
will return to it shortly, after gaining some confidence in the power of the 
matrix model. 
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7.1 Comparison for b~ t = = b 

Under this assumption the last term in equation (jll2j) simplifies considerably, 
and the effective superpotential reduces to: 

V , S\ „. detM (6m6)detM detm . 

W eff = S 1 - In — - 5 In - 9(N - + 5 In— — . 113 

V A 3 / A 2 ^- 1 ) detm A 2 

To compare with the gauge theory effective superpotential we must integrate 
out the glueball superfield. 

^f = -^ = ^ (114) 

5S A detm detM 

and thus the effective superpotential is given by: 

A 2iVc_1 detm (bmb) , - 

W eff = s ^ i det M , 115 

detM detm 

which reproduces the field theory result. The first term can be recognized as 
the ADS superpotential upon noticing that A 2Nc ~ 1 det m is the scale of the 
theory obtained from the initial one by integrating out two quarks with mass 
matrix m. 



7.2 General analysis 



We now turn to analyzing the last term in equation ()112|) and isolating the 
part arising from planar and single boundary (in the sense of j22]) diagrams. 

It is easy to reorganize this term using equation (|l()9jl . To avoid cluttering 
the equations, let us introduce 



A = detm - (bM~ l b) det M 
Then the last term in ()109|) becomes: 



;ii6) 



In det (5 b a det m + b t PV^V 3 



cb 



TrmA5° + Trln 



K + jX1Q° 3 QI 



S In A + In det 



(117) 



where, as before, we identified the 't Hooft coupling with the glueball su- 
perfield and the matrix whose determinant is computed in the second term 
carries flavor indices. 
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We must now identify the leading terms in this equation - terms generated 
by planar diagrams with as many gauge index loops as the diagrams with 
one boundary. For this purpose it is important to notice that the compu- 
tations in the previous section yield the sum of all 1-loop n-point functions. 
Furthermore, the planar, 1-boundary contribution must be proportional to 
the number of colors iV c = S, since there is one gauge index loop in such dia- 
grams. It is therefore clear that only the first term in equation (|117|) should 
be kept since the determinant in the second term is in flavor space and there 
is no term in its expansion which is proportional to the number of colors. 
Thus, the gauge theory effective superpotential is given by: 



W. 



eff 



Sin 



detM 

A2(Jv e -i) 



118) 



(bmb)detM 1 /, „ « . 

i k — = ^ + S In — ( det m - (6M _1 6) det M 

det m - {bM-^b) detM A 2 V 



Integrating out S leads to: 



S_ 

A3 



A 2(7V c -2) 



det m 



(bM' l b) 



.detM 

which in turns implies that the effective superpotential is: 



119) 



A 



2N c -l 



det m 
detM 



- (bM~ L b) 



(brab) det M 



detm- (6M- 1 6) detM 



(120) 



This reproduces the field theory result (j74[ 198 )1 . 



8 Vacua 

Although it is already clear that there is an exact agreement between the 
matrix model and gauge theory, let us briefly discuss the vacua of the gauge 
theory and their construction from the matrix model. In gauge theory we 
need to integrate out all mesons and baryons, while on the matrix model 
side we need to compute the full partition function. We begin with the 
gauge theory discussion. We will discuss the construction in the language of 
section El and relate it at the end with section El 
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8.1 Integrating out all fields in gauge theory 

Let us recall equation (|TE|). which determines the low energy effective su- 
perpotential for a general Nf 7^ 3: 

(N f - lfx-y^-^xf{x)[f{x) + (N f - l)(N f - 3)xf'{x)] = -1. 

Besides this equation, there are other equations f(x) obeys, obtained from 
varying the dynamical superpotential with respect to the mesons: 

B'Bj - (detM)(M- 1 )' i i 

Using various equations from section 3.1 this equation can written as: 
- (det M)(M" 1 )} = -m)k 2Na - x 

= ^^k (N f - l) 2 \(N f - l)V / (x)[(detm)A 2 ^~ 1 ] 2 /W- 1 )^S 
(bmb) L (bmb) 

+A 2N ^(f(x) - 2(N f - l)xf(x)) N f- 2 (f(x) + (N f - l)(N f - 3)xf'(x)) 

- m )\ 2N <- l (f(x) - 2(N f - l)xf'(x)ff- 2 [f(x) + (N f - l)(iV> - 3)xf'(x)} . 

(122) 

To satisfy this equation, the coefficient of (mb)i(bm)j in the square bracket 
must vanish, and the coefficient of m* must agree on both sides. Therefore 
we find 

(N f - lYx-Wt-Vxfix) + (f(x) - 2(N f - l)xf'(x)) N f- 2 = ,(123) 
(f(x) - 2(N f - l)xf'(x)) N f- 2 (f(x) + (N f - l)(N f - 3)xf'(x)) = 1 .(124) 

Thus, there seem to be three equations for a single function; it turns out 
however that one of them can be obtained from the other two. In general, 
we cannot expect to find a consistent solution for two first-order differential 
equations for one function. In the Nf = 3 case, the two differential equa- 
tions were self-consistent, and their combined effect was to fix the integration 
constant in f(x) 9 . We expect the same to happen here. 

In general, we can solve for f'(x) using Eq. (JTHJ), and substitute it to one of 
the other equations. Since Eq. (|T6|) is quadratic in f'(x), it has two solutions. 

9 Indeed, if one did not fix the integration constant c in f(x) = ±V c — x 1 / 2 using the 
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Only one of them is consistent with the boundary condition |/(0)| = 1. 
Keeping only the consistent solution, we find 



'(Nf - l)(jV) - 2)f - y/(N f - Iff 2 - A{N f - l)(N f - 3)z 
(N f -3)(N f -l) 



z- 1 (-(N f - 1)/ + y/(N f - lYP - 4(iV / - l)(N f ^3)~z) 

2(N f - 3) 



N f ~2 



0,(126) 



where z = x 1 '^* -1 '. This equation determines the function f(x) implicitly. 

The same results can be obtained following the steps in section El In 
particular, when all Nf = N c + 1 flavors are integrated out, equations (J93|) 
and (jUl)) become: 

detM = -BX-— (127) 
B 

detM^- 1 +XBdetM N f~ 2 = det m A N f( 2N ^ (128) 

The effective superpotential is then obtained by substituting the solutions of 
these equations in the superpotential (jHOjl 

A 2Nc ~ 1 W e ff = (B 2 X + N f - 1) det M + B 3 X 2 + 3BX. (129) 

It is not hard to check that this reproduces the results in chapter 3 for the 
case of an SU(2) gauge group with 3 flavors; it is, however, somewhat more 
challenging to see that it agrees with ()126p as well. 

8.2 The Matrix Model Free Energy 

Let us now consider the matrix integral we considered before, but with all 
flavors massive. In this case, we can reinterpret the 5-function as arising from 
the change of variables 

/ DQDQ = J DM J DQDQ5(QQ - M) (130) 



boundary conditions, equation H124|l : 



l T -l/2 



2VC-X 1 / 2 I 



would fix this constant to be c = 1. 
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Thus, to find the effective superpotential as a function of the glueball su- 
perfield (J2J), we must supplement the results of the previous section with a 
mass term for the remaining mesons and then compute the integral over M 
as well. We recall that we are interested only in the 1-boundary free energy. 
Thus, the integral can be computed by a saddle-point approximation. Al- 
ternatively, it is easy to see that the one-boundary free energy is given by 
the sum of all tree-level Feynman diagrams arising from the superpotential 
()118|) . This implies that, as expected, the effective superpotential is unique 
even when expressed in terms of the glueball superfield. The vacua of the 
theory arise in this language as the critical points of W e ff(S). 

We now illustrate this simple observation for the SU(2) theory with three 
flavors, leaving to the reader the exercise of recovering the more involved 
results of section IHTT1 

8.3 Back to 5*7(2) 

Consider the equation (|118|) for the case of an 577(2) theory with three flavors. 
Since M is 1-dimensional, the superpotential is: 

W eff = 5(l-ln^)-51n^ (131) 



(bmb)M 



+ S In — (det m - (6,6,)) + mM 



det m - (6,6,) ' A 2 

where b\ and b\ are the sources with indices along the meson which was not 
integrated out in the previous section. The saddle point equation is: 

S (bmb) det m- (bmb) v /m^ 0\ 

— = m — = = — where m — „ 132 

M detm-(6,6,) detm-(6,6,) \0 m) 

and b and b are understood as 3-component vectors. Then, the effective 
superpotential as a function of the glueball superfield is: 

W eff (S) = 25(l-ln^)-51n A ' (133) 
V A J / det m — {bmb) 

As argued before, the vacua are now described by the critical points of 
Weff(S), and are given by 

S , det m — (bmb) „ / ; =~ , , /9 

21n— = ln '- S = ±y/detm- (bmb) A 3/2 (134) 
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The superpotential at the critical points is therefore: 

W eff \ = ±2 V / detm- (bmb)A 3/2 . (135) 

We thus recover the matrix model result found directly in equation (J27|) of 
section |U as well as the field theory result. 



9 Baryons and Geometric Transitions 

In this section we discuss the baryons in the context of the geometric tran- 
sitions. The gauge theory is engineered by wrapping D5 branes on several 
compact P 1 cycles of a geometry which locally, around each cycle, is the 
geometry of the small resolution of the conifold. Alternatively, it can be de- 
scribed using the T-dual brane configuration, where the D5 branes wrapped 
on P 1 cycles are mapped into D4 branes stretched between NS branes [TH] . 

I5HG21G3]. 

Let us begin by briefly reviewing the results of jHU EH] , concerning the 
baryonic degrees of freedom in MQCD. First, we need to comment on having 
an SU(N) rather than an U(N) gauge group. The Type IIA brane config- 
uration as well as the Type IIB geometric construction describe a classical 
U(N) gauge theory. The M theory limit describes a quantum SU(N), where 
the U(l) factor decouples. As explained in [5T] . the U(l) factor is recovered 
after the geometric transition, when the SU(N) part confines. Therefore, the 
approach of [TH] cannot be applied for the case of baryons, as the quantities 
in matrix models were obtained from the parameters of brane configurations 
via lifting to M theory. 

It is nevertheless possible to collect some information about the vacuum 
expectation values of the baryon operators in MQCD. As described in 
in the case Nf = N c , the difference between a baryonic and a non-baryonic 
branch is that the asymptotic regions of the former intersect, and the ones 
of the latter do not. Indeed, the asymptotic regions for the non-baryonic 
branch are given by: 

t = (w 2 + A% =l ) N °' 2 , v = 

t = A^ c ! , w = (136) 

while the ones for the baryonic branch are given by: 

* = w 2Nc , v = 

t = A 2 Jti , w = . (137) 
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It is clear that the two branches intersect in (|136|) , but are separated in (|137|) . 
The distance between the asymptotic regions in (|137|) is the value of BB. 

In M theory terms the geometric transition corresponds to a transition 
from an M5 brane with a worldvolume containing a Riemann surface in the 
(v,w,t) plane to an M5 brane with two dimensions embedded in (v,w), for 
constant t. The equation in (v , w) represents an NS brane which is T-dual 
to the deformed conifold. 

In the case of fJ136|) - ()i37j) . v and w are decoupled so the above discussion 
does not apply. In the language of [THj, this can be understood by starting 
with D4 m branes corresponding to massive flavors, taking the mass to zero 
and combining with a color D4 brane to get a D4^ brane which describes 
a flavor with an expectation value. Therefore, in the geometrical picture, 
there are no D5 branes on the compact P cycles and there are only D5 
branes on the noncompact 2-cycles. We then see that the duality between 
matrix models and field theory fails in this case. 

The only way to use the results of EJ |3] is to give mass to one of 
the flavors, which means decomposing one D4jv/ brane into a D4 m brane 
and a color brane. This is exactly the procedure discussed in detail in [32] 
where a method to deal with this case was stated. Therefore, we see that 
the difficulties with the matrix model analysis of the baryon operators have a 
geometric counterpart. This should probably be expected, since the geometry 
is underlying the matrix models. 

10 Conclusions 

In this paper we further analyzed the extension of the Dijkgraaf-Vafa proposal 
to theories containing fields in the fundamental representation. While this 
extension was thoroughly analyzed in situations in which the gauge theory 
was described solely in terms of mesons, the matrix model description of 
baryonic deformation remained until now largely unexplored. The main goal 
of our work was to fill this gap. 

We have started with the Af = 1 SQCD with gauge group SU(N C ) and 
N c + 1 flavors whose effective superpotential was conjectured in [2] and de- 
formed the theory by adding baryon sources as well as mass terms for either 
two or all flavor fields. We compared the resulting effective superpoten- 
tial obtained by integrating out the appropriate mesons and baryons with 
the one coming from the matrix model computations and we found perfect 
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agreement. Of essential importance has been the correct identification of 
Feynman diagrams contributing to the superpotential. 

We expect that the effective superpotential for other s-confining theories 
is computable using matrix model techniques along the lines described here, 
after suitable deformations by mass terms and other sources. 

SQCD theories with Nf > N c + 2 are usually analyzed using Seiberg's 
duality. One may ask whether the matrix model techniques can shed light 
on their effective superpotential. Using 't Hooft's anomaly matching condi- 
tions it was shown that the mesons and baryons are not the only low energy 
degrees of freedom. However, the complete set of low energy fields is not 
known. Nevertheless, by inserting sources for the known fields in the tree 
level superpotential, the matrix model perturbation theory should allow one 
to recover the truncation of the full effective superpotential to these fields. 
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